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EXECUTIVE  OVERVIEW 


Many  attributes  of  port  use  are  random  in  nature  and,  as  such,  complicate 
the  problems  associated  with  port  planning  and  design.  This  work  recognizes 
this  situation  and  describes  a  method  by  which  statistical  techniques  can  be 
used  to  analyze  port  resource  requirements.  The  methodology  derived  is  titled 
Port  Systems  Requirements  Prediction  Methodology.  Its  description  is  intended 
to  provide  the  information  required  by  the  full  range  of  potential  users  -  from 
manager  to  technician.  A  sample  application  is  described  as  a  mechanism  to 
display  data  file  development  and  computer  programs. 

The  premise  of  the  technical  approach  is  that  one  cannot  look  into  the 
future  for  design  information,  but  rather  the  statistical  distribution  of 
utility  service  requirements  can  be  established.  These  distributions  can  then 
be  used  to  size  port  systems  so  that  a  predetermined  low  probability  of  failing 
to  meet  a  demand  can  be  achieved. 

However,  the  Port  Systems  Requirements  Prediction  Methodology  is  not 
restricted  to  a  single  analytical  output.  This  report  discusses  its  potential 
use  as  a  means  of  communication  for  MILCON  planning  and  as  a  vehicle  for 
analyzing  the  output  of  port  utility  measurement  projects.  In  short,  the 
Port  Systems  Requirements  Prediction  Methodology  is  a  powerful  analytical 
tool  and  the  challenge  presented  by  it  is  for  planners  and  designers  to 
articulate  problems  and  issues  that  can  exploit  its  potential. 
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SECTION  ONE 
INTRODUCTION 


1 .0  Background. 

The  Port  System  Project  Program  Plan  established  the  objective  of 
relating  actual  utility  demand  data  to  pier  design  criteria.  There  are  two 
reasons  for  doing  this.  The  first  is  to  update  the  criteria  used  to  size  the 
utility  systems  on  Navy  piers.  And  secondly,  it  is  intended  to  be  able  to 
predict  pier  system  changes  that  will  be  requirr  by  changing  fleet  demands. 

A  prediction  technique  such  as  this  would  be  a  valuable  tool  for  port  master 
planners  to  use  when  evaluating  various  design  options.  The  work  described  in 
this  report  produced  a  model  that  relates  the  demand  for  services  at  a  selected 
port  facility  to  the  operating  tempo  of  the  ships  assigned  to  that  facility  and 
individual  ship's  cold-iron  requirements.  The  model  includes  computer  simula¬ 
tions  with  enough  iterations  for  statistically  valid  conclusions. 

1.1  Objective. 

The  objective  of  this  work  is  to  develop  a  methodology  for  predicting 
port  service  requirements  for  a  given  port  load. 

1 .2  Scope. 

The  methodology  described  in  this  report  covers  the  concept  formula¬ 
tion,  the  development  of  representative  data  files,  computer  programming, 
simulation  and  illustrations  of  potential  interpretations  of  results. 

1 .3  Approach. 

The  approach  to  developing  the  Port  System  Requirements  Prediction 
Methodology  starts  with  the  premise  that  the  need  for  port  services 
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can  be  analyzed  by  calculating  the  utility  demands  of  specified  ship  popula¬ 
tions  while  recognizing  the  service  demands  of  individual  ship  classes.  This 
information  is  combined  with  representati ve  ship  movement  data  to  form  a  port 
activity  model  that  can  be  computer  simulated.  This  in  turn  permits  the  proba¬ 
bility  distribution  of  utility  demands  for  specified  ship  populations  to  be 
statistically  estimated.  It  is  this  information  then,  the  probability  distri¬ 
bution  of  the  utility  demands  of  specified  ship  populations,  that  constitutes 
the  primary  output  of  the  Port  Systems  requirement  methodology.  Formats  for 
aggregating  and  interpreting  output  information  are  suggested  as  the  final 
step  in  methodology  formation. 

1 . 4  Limi tations. 

The  work  described  herein  represents  the  development  of  a  methodology 
only.  The  results  are  based  on  representative  data  and  care  must  be  exercised 
not  to  use  the  results  of  this  study  out  of  context.  For  use  of  the  methodo¬ 
logy  beyond  the  demonstration  described,  the  data  files  would  require  adjust¬ 
ment  for  the  specific  situation.  In  addition  to  having  used  representative 
data,  the  following  simplifying  assumptions  have  been  made: 

-  The  only  utility  analyzed  was  electric  power,  and  each  ship  was 
assumed  to  be  drawing  its  full  design  capacity  for  shore  power 
in  accordance  with  DM-25. 

-  No  consideration  was  given  to  the  time  of  day  variations  that  would 
be  expected  for  electric  power  demand. 

-  The  number  of  ships  at  cold-iron  in  each  ship  class  is  assumed  to 
be  binomially  distributed  and  independent  of  other  classes. 

-  Resource  demands  per  ship  do  not  correlate  with  the  number  of  ships 
present. 

It  should  be  noted  that  these  assumptions  are  not  a  prerequisite  to  the  correct¬ 
ness  of  the  technical  approach  to  the  methodology  development.  They  have  been 
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made  to  limit  the  scope  of  this  work  and  to  provide  a  reasonable  example  of 
its  application.  In  fact,  each  limitation  can  be  eliminated  by  the  use  of 
more  precise  information.  Methodology  products  would  then  provide  an 
increasingly  accurate  representation  of  waterfront  activity. 

1.5  Application. 

The  immediate  application  of  this  methodology  is  to  establish  the 
cummulative  probability  distribution  for  utility  demand  for  homeport  loads  of 
varying  size  and  ship  class  mixes.  This  would  be  useful  for  sizing  utility 
system  improvements  as  well  as  for  estimating  the  time  to  initiate  such 
projects.  Another  application  would  be  to  use  the  Port  System  Requirements 
Prediction  Methodology  to  interpret  the  output  of  port  utility  sampling  projects. 
In  this  way,  the  results  of  this  study  can  serve  the  Port  Systems  Project  by 
providing  a  focus  for  the  various  utility  demand  data  experiments  being  planned. 
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SECTION  TWO 


METHODOLOGY  CONCEPT 


2.0  Objecti ve 

The  objective  of  this  section  is  to  show  how  to  use  a  statistical  approach 
for  port  requirements  planning  and  to  introduce  a  method  by  which  this  may  be 
done.  This  method  is  called  the  Port  System  Requirements  Prediction  Methodology. 

2.1  Scope 

This  section  discusses  the  broad  qualities  of  the  port  planning  problem, 
describes  a  methodology  for  analyzing  typical  waterfront  problems  using  statis¬ 
tical  techniques,  and  presents  potential  uses  for  the  products  generated  by 
the  Port  Systems  Requirements  Prediction  Methodology. 

2.2  Discussion 

Ship  service  requirements  at  a  given  port  can  be  predicted  accurately  for 
only  a  few  days  or  weeks.  This  is  due  primarily  to  the  uncertainty  in  future 
fleet  operating  characteristics  which  directly  impact  the  presence  or  absence 
of  any  given  ship  at  the  port  at  any  specified  future  time.  There  are  also 
uncertainties  in  the  per  ship  demand  and  in  the  number  and  type  of  ships  opera¬ 
ting  from  a  given  port.  The  result  is  that  service  capacity  decisions  must  be 
made  with  limited  information  regarding  operational  scenarios.  It  cannot  be 
stated,  for  instance,  that  on  a  specific  date  in  1990  there  will  be  83  ships  in 
San  Diego  Harbor.  Hence,  maximum  resource  requirements  cannot  be  established 
by  surveying  a  list  of  ships. 
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Under  these  ci rcumstances  a  reasonable  approach  for  port  system  analysis  is 
to  attempt  to  derive  the  probability  distribution  for  resource  demand  rates  for 
different  homeporting  plans.  The  capacity  for  a  selected  utility  can  then  be 
represented  as  a  preset  low  probability  of  failing  to  meet  the  demand  during  the 
design  lifetime  of  the  capital  installation.  This  approach  offers  reasonable 
possibilities  for  validation  by  reviewing  historical  data  and/or  through  planned 
observations  over  the  early  years  of  its  application. 

Figure  2-1  illustrates  how  information  of  this  nature  could  be  presented: 


DEMAND 

FIGURE  2-1  DESIRED  FORMAT  FOR  METHODOLOGY  OUTPUT. 
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The  abscissa  is  the  demand  level  for  a  given  resource.  The  ordinate  is 
the  probability  that  the  requirement  of  a  specified  group  of  ships  for  that 
resource  will  exceed  a  specified  demand  level  during  an  interval  of  interest. 

The  interval  can  be  defined  in  many  ways;  it  may  be  a  randomly  selected  day  in 
a  month,  part  of  a  day,  an  entire  year,  or  even  the  lifetime  of  the  capital 
investment.  Figure  2-1  depicts  two  graphs  for  the  same  interval.  Each  curve 
represents  a  different  port  load.  One  might  be  for  the  current  homeport  plan, 
and  the  other  for  some  future  port  load.  In  this  case,  the  preset  "failure" 
probability  Pj-  is  established  relative  to  the  present  supply  level,  S.  When 
the  demand  reaches  the  supply  level  S,  the  system  will  no  longer  function  without 
special  management  measures  such  as  rationing,  berthing  shifts,  lighting  off 
ships  utility  systems,  etc. 

The  relationship  between  the  curves  of  Figure  2-1  enables  one  to  quickly 
evaluate  alternative  strategies  for  dealing  with  the  resource  requirement.  For 
example,  if  there  is  no  change  in  the  supply  capacity,  then  the  probability 
of  failure  would  rise  to  the  value  Pp'  in  the  case  of  the  set  of  ships  in  the 
future.  Alternatively,  if  it  is  required  that  the  failure  rate  remain  constant, 
then  the  supply  would  have  to  increase  from  S  to  S  .  A  larger  increase  in 
supply  would  result  in  a  decreasing  failure  probability  while  either  a  zero  or 
a  small  increase  would  result  in  a  larger  probability  of  failure.  This  infor¬ 
mation  can  be  of  immediate  use  in  formulating  priorities  for  acquisition  of 
waterfront  improvements. 

2.3  Port  System  Requirements  Prediction  Methodology  Objective 

Developing  information  similar  to  that  depicted  by  Figure  2-1  for  various 
resources  constitutes  a  new  technique  for  analyzing  port  requirements  and  is  the 
objective  of  the  Port  System  Requirement  Prediction  Methodology.  This  technique 
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if  valuable  from  several  points  of  view.  For  example,  the  degree  to  which 
a  port  control  authority  is  forced  to  exercise  special  management  procedures 
because  of  utility  shortages  can  be  quantitatively  presented  to  Fleet  decision 
makers.  Alternatively,  the  impact  of  not  updating  pierside  resources  to 
accommodate  planned  ship  alterations  or  new  ship  acquisitions  can  be  quickly 
determined.  Overall,  the  information  available  through  this  methodology  is 
obtained  at  low  cost  and  in  a  short  timeframe  once  the  methodology  has  been 
finalized.  Other  advantages  in  terms  of  project  justification  and  acquisition 
strategies  can  be  developed  if  this  methodology  is  placed  in  the  hands  of 
knowledgeable  port  planners. 

2.4  Methodology  Application 

The  probability  distribution  of  resource  demand  rates  (i.e.,  the  amount  of 
the  resources  that  would  be  consumed,  if  available)  is  identified  above  as  in¬ 
formation  that  could  be  effectively  used  for  port  requirements  analysis.  Graphs 
of  the  nature  of  Figure  2-1  are,  then,  the  goal  product  of  the  Port  Systems 
Requirement  Methodology.  Hence,  the  procedure  for  developing  this  information 
is  of  primary  interest,  and  is  in  fact  the  conceptual  basis  of  the  Port  System 
Requirements  Methodology.  Figure  2-2  is  a  model  of  how  the  goal  product  of 
Port  System  Requirements  Methodology  is  generated.  This  figure  shows  that  there 
are  three  phases  of  the  Methodology.  First  is  data  acquisition  and  assembly. 

In  the  example  used  to  illustrate  the  methodology  6QHz  power  was  selected  as  the 
resource  to  be  analyzed.  Ship  demand  data  were  taken  from  DM-25.  (These  data 
could  be  made  much  more  realistic  by  using  actual  measurement  data  from  the 
ports  of  interest).  Ship  populations  and  cold  iron  probabilities  were  derived 
from  OPNAV  INST  3111. 14U  and  NAVSEC  Report  6139-72-2(1977),  respectively. 
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OPNAV  f NS r  512 1 . :-'J 

NAVSEC  PEPT  6139-72-7  <1977) 
PS  MASTER  PLAN 
SEWELLS  POINT 


PSP. FOR  IS  THE  NAME  GIVEN 
TO  THE  PROGRAM  THAT  COMPUTES 
THE  DEMAND  DISTRIBUTIONS 
BY  BOTH  THE  CENTRAL 
LIMIT  THEORM  AND  BY 
SIMULATION.  A  COPY  OF 
PSP. FOR  IS  IN  APPENDIX  A. 


PLTF11  AND  PLTFI2  ARE  THE 
PROGRAMS  PREPARED  TO  PROVIDE 
AUTOMATIC  FORMATTING.  COPIES 
ARE  IN  APPENDIX  A. 


FIGURE  2-2.  ACTIVITY  FLOW  IN  PORT  SYSTEMS  REQUIREMENTS  PREDICTION  METHODOLOGY 
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Next,  Naval  Station,  San  Diego  and  Naval  Operating  Base,  Norfolk  were 
selected  as  the  ports  under  study.  The  second  phase  consisted  of  the  actual 
generation  of  the  goal  product,  i.e.,  the  demand  distribution.  Figure  2-2 
shows  that  two  techniques  were  used  to  develop  this  produc£,  the  computer 
simulation  and  the  Central  Limit  Theorem  methods.  Computer  programs  were 
prepared  during  this  study  which  performed  the  required  calculations  and 
permitted  automatic  graphical  presentation  of  the  results.  The  last  phase 
of  this  methodology  was  output  analysis.  This,  and  the  mathematical  theory 
that  forms  the  basis  of  computational  work  are  described  in  the  following 
sections  of  this  report. 
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THEORETICAL  BASIS 


3.0  Objective. 

The  objective  of  this  section  is  to  describe  the  theoretical  basis  of 
the  Port  System  Requirements  Prediction  Methodology. 

3.1  Scope. 

Three  distinct  elements  form  the  theoretical  basis  of  the  Port  System 
Requirements  Prediction  Methodology.  These  are: 

a.  An  interpretation  of  port  conditions  that  can  be  used  as  a 
framework  for  statistical  analysis. 

b.  The  derivation  of  the  mathematical  expressions  used  to  generate 
the  probability  distribution  for  utility  demand. 

c.  An  explanation  of  the  alternate  means  of  acquiring  the  products 
of  the  Port  System  Requirements  Prediction  Methodology. 

The  first  two  elements  are  discussed  in  a  single  subsection  titled 
"Theoretical  Approach"  and  the  latter  element  is  discussed  under  the  subsection 
titled  "Product  Development." 

3.2  Theoretical  Approach. 

The  distribution  of  total  demand  on  a  given  service  facility  will 
depend  on  a  number  of  port  conditions.  The  most  significant  of  these  conditions 
are: 

1.  The  ship  population  that  may  use  the  service  facility. 

2.  The  probability  that  a  member  of  each  individual  ship  class  will 
actually  be  at  cold  iron,  and 
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3.  The  probability  distribution  of  the  demand  for  each  individual  ship 
given  that  is  at  cold  iron.  In  general  terms,  this  may  be  a  joint 
probability  distribution  for  more  than  one  ship  if  there  is  a 
cooperation  of  some  type  between  ships. 

The  use  of  these  conditions  to  determine  the  distribution  of  the  total  demand 
for  a  specific  utility  is  described  below. 

First,  let  Pk  be  the  probability  that  a  ship  in  the  kth  class  is  not  at 
cold  iron  in  the  port  or  service  facility  of  interest,  and  let  Nk  be  the  total 
number  of  ships  in  class  k  that  are  assigned  to  the  port  under  study  during 
the  period  of  interest.  The  cold  iron  probability,  P^,  is  simply 


P  ‘  P 
iK  y2k 


where 


=  the  probability  that  a  ship  in  the  kth  class  will  be 
present  at  the  port,  and 

=  tiie  probability  that  a  ship  in  port  in  the  kth  class 
will  be  at  cold  iron. 


The  next  step  is  to  estimate  how  the  numbers  of  ships  present  are  statis¬ 
tically  distributed.  As  an  initial  step,  we  assume  that  there  is  no  correlation 
between  the  presence  of  any  one  ship  and  any  other  ship.  This  will  not  be  true 
if  ships  are  traveling  in  company,  but  is  a  reasonable  initial  representa¬ 
tion.  In  this  event,  the  number  of  ships,  n^,  actually  present  in  the  interval 
of  interest  will  be  binomially  distributed.  The  related  probability  is: 


P<V  ■ 


Mp  (l  -  pJ  Nk'"k 
V  K  k 


where  the  first  term  on  the  right  hand  side  is  the  number  of  combinations  of  Nk 
articles  that  can  be  taken  nk  at  a  time,  i.e. 

/NkU  _ Uhl _  (2) 

Vnk /  nk:  (\  -  nk); 
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Figure  3.1  presents  an  example  of  a  binominal  distribution.  The  total  population, 
Nk ,  has  been  set  at  5  and  the  cold  iron  probability  has  been  set  to  0.5.  The 
probability  distribution  obtained  for  these  parameters  is  the  familiar  coin¬ 
toss  distribution  giving  the  probability  of  obtaining  a  given  number  of  heads 
(or  tails)  when  5  coins  are  tossed.  It  should  be  noted  that  although,  in  the 
example  given,  the  distribution  is  symmetric,  this  will  not  usually  be  the  case. 


The  mean  or  average  value  of  n^  is  denoted  by  E  (n^)  indicating  "expected 
value"-  It  is  relatively  easy  to  show  that  the  mean  has  the  following  values: 


E<\>  ■  \pk 


(3) 


Additionally,  we  define  the  variance  V  (n  )  of  the  distribution  relating  to  the 
dispersion  about  the  mean.  It  is,  in  fact,  equal  to  the  average  of  the  squares 
of  the  deviations  from  the  mean  value.  For  the  binominal  distribution,  the 
variance  is  given  by: 


v‘"k'  ■  W'-V 


(4) 


In  the  example  given  above  the  expected  va  lue  is  5  x  0.5  =  2.  5  (ships)  while  the 

2 

variance  is  1.25  ships  .  The  standard  deviation  of  the  number  of  ships  present 
is  defined  as  the  square  root  of  the  variance  i.e.,  \j  1.25  or  1.12. 

The  demand  for  resource,  j,  by  a  given  ship  will  depend  on  its  class,  k. 

If  dk  is  the  value  of  this  demand,  it  will  be  distributed  over  a  range  depending 
on  the  operations  undertaken  during  servicing.  Thus  for  a  given  interval  of 
time  there  will  be  a  matrix  of  probabilities  Pikj  giving  the  probability  that 
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Total  population  in  ship  class  of  interest 


p  =  Probability  that  class  member  is  at  cold  iron 


Probability  that  n^  ships  are  present 


Mean  number  of  ships  present  =  N^p 
Variance  in  number  of  ships  present  =  N^p(l-p) 

RGURE  3-1.  Example  of  a  Binomial  Distribution 
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djk  will  be  in  the  range,  d j k ( i - 1 )  t0  djk(i)  »  where  djk(i)  is  the  marker 
for  the  given  resource.  In  testing  our  program,  however,  we  have  made  the  assump¬ 
tion  that  variability  in  the  per  ship  demand  is  less  important  than  the  varia¬ 
bility  in  the  number  of  ships  present  and  have  taken  the  value  of  each  d  to 

JK 

be  a  constant.  This  is  equivalent  to  having  a  single  infinitesimal  range  and 
setting  P  to  unity  for  all  j  and  k.  With  this  assumption,  the  mean  or 
expected  total  demand  for  a  given  resource  and  its  variance  can  be  estimated 
as  follows: 
i s  gi ven  by 


as  follows:  If  D.  is  the  total  demand  for  resource  j,  then  its  average  value 


E(D  )=  EV  E(nk)  •  EWjk 
k  k 


(5) 


and  the  variance  is  given  by 


-  XXjk  «\>  -  &jk  Vk  u-pk)  (6) 

k  k 

It  is  equations  (5)  and  (6)  then,  that  are  used  to  generate  the  goal  product 
of  the  Port  Systems  Requirements  Prediction  Methodology. 

3.3  Product  Development 

The  calculation  of  the  distribution  of  the  total  demand  D.  for  resource  j 
is  quite  complicated  and  is  best  performed  by  multiple  complementary  approaches. 
Fortunately,  the  three  methods  that  have  been  developed  always  provide  a 
mechanism  for  applying  two  techniques  to  every  case  considered.  We  shall  now 
present  the  details  of  the  three  methods. 


Method  1  -  Explicit  Enumeration 

This  approach  is  suitable  for  "small"  problems,  i.e.  when  the  number 
of  ships  and  ship  classes  are  small.  The  probability  of  obtaining  a  given  dis¬ 
tribution,  ,  of  ships  in  class  k  is  given  by 

m 

P(r\y  n2,  ...nm)  =  P(n1)P(n2) . .  .P(nm)  s  IT  p(nk)  (7) 

k  —  1 
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In  other  words,  the  probability  that  we  would  obtain  n^  ships  in  class  1, 
n^  ships  in  class  2,  etc.,  is  simply  the  product  of  the  probabilities  that  the 
states  be  obtained  independent  of  the  other  classes.  If  we  believe  that  the 
numbers  of  members  of  any  class  do  not  depend  on  the  numbers  present  in  any 
other  classes,  but  only  on  the  total  population  in  that  class  and  the  cold 
iron  probabilities,  then  the  use  of  the  product  with  be  valid  with  the  indivi¬ 
dual  probabilities  given  by  equation  (1). 

The  number  of  terms  in  equation  (7)  rises  rapidly  as  the  harbor  size 

grows.  If  Nk  is  the  number  of  ships  in  class  k,  there  are  (Nk+1)  possibilities 

for  this  case.  One  is  added  to  allow  for  the  possibility  of  zero.  The  total 

m 

number  of  terms  is  n  (N.+l).  As  an  example,  if  there  were  ten  classes,  each 
k=l  K 

with  ten  ships,  then  the  number  of  possible  populations  in  the  port  would  be 
ll10.  Thus  while  this  approach  is  useful  for  validating  calculations  on  small 
problems,  it  is  normally  impracticable  for  realistic  cases  and  has  not  been 
formally  programmed. 

Method  2  -  Monte  Carlo  Simulation 

Here,  we  introduce  the  concept  of  a  "realization"  of  the  conditions  of 
operation  of  the  port.  If  we  are  given  the  data  describing  the  probability 
that  the  port  will  be  in  any  given  condition,  then  through  the  use  of  a  random 
number  generator  we  can  perform  a  computer  operation  simulating  that  condition 
with  a  likelihood  of  occurrence  equal  to  the  real  world  probability.  This  is 
a  realization  of  this  condition.  By  generating  many  such  realizations  we  can 
compute  the  statistics  of  usage  of  any  resource  surveyed. 
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For  each  class  we  establish  a  random  number  of  ships  at  cold  iron  in 
accordance  with  a  binomial  distribution.  Figure  3-2  illustrates  how  this  is 
done.  We  generate  a  pseudorandom  number,  p,  rectangularly  distributed  in  [0,1] 
and  then  identify  the  number  of  ships  present  to  the  largest  number,  n,  such 
that  the  probability  of  obtaining  fewer  than  n  numbers  is  less  than  p.  For 
purposes  of  this  work,  a  pseudorandom  number  is  considered  to  be  a  number  in 
the  range  [0,1]  that  is  generated  by  a  computer  program  and  whose  occurrence 
follows  a  rectangular  distributon  so  as  to  possess  the  statistical  quality  of 
randomness.  For  Figure  3-2,  the  magnitudes  of  the  probability  increments  are 
equal  to  the  probabilities  given  in  Figure  3-1  and  in  this  specific  example, 
the  number  selected  by  the  procedure  given  is  2. 

The  generation  of  a  realization  of  the  entire  system  involves  the  deter¬ 
mination  of  the  nuirber  of  ships  present  in  each  class.  This  requires  that 
the  operation  described  above  be  performed  once  for  each  class.  Having  deter¬ 
mined  the  number  of  ships  present  in  each  class,  we  can  then  determine  the 

total  usage  of  resources  for  this  realization. 

% 

The  goal  product  is  accomplished  by  generating  a  large  number  of  realiza¬ 
tions  of  the  entire  system  (at  least  100)  and  counting  the  number  of  times  that 
the  total  demand  for  any  given  resource  exceeds  a  specified  value.  An  advantage 
of  this  approach  is  that  complications  such  as  a  probability  distribution  of 
resource  demands  by  individual  ships  can  readily  be  accommodated  in  this  part 
of  the  calculation. 
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NO.  OF  SHIPS,  n 


•  RECTANGULARLY  DISTRIBUTED  PSfc'UDO— RANDOM  NUMBER 
<0<A< 


FIGURE  >2.  METHOD  OF  CALCULATING  NUMBER  OF  SHIPS  AT 
COLD  IRON. 


It  must  be  remembered  that  the  accuracy  of  the  simulation  result  improves 
the  number  of  realizations,  N,  is  increased.  In  general,  the  standard  error 
p,  in  any  probability  estimate,  p,  is  given  by  \N  p(l-p)  so  that,  in  order 
to  halve  the  error,  we  must  increase  the  number  of  simulations  by  a  factor  of  4. 

Initially,  we  used  a  value  of  N  of  100,  but  this  number  can  be  increased 
as  the  input  data  becomes  more  accurate. 

Method  3  -  Central  Limit  Theorem 

The  total  demand  for  any  given  resource  often  arises  from  many 
randomly  behaving  components  that  have  the  same  probability  density  functions, 
and  is  computed  as  a  sum  of  these  components.  The  Central  Limit  Theorem  states 
that,  under  these  conditions,  the  sum  is  distributed  in  an  asymptotically 
"Normal"  or  "Gaussian"  fashion  independent  of  the  probability  distributions  of 
the  individual  components. 

The  normal  distribution  for  any  variable,  x,with  mean  and  variance  E(x) 
and  V(x),  respectively,  takes  the  form: 

P(x)  =  (V2*v<*f  )*’  exp  .  (8) 

and  the  integrals  of  this  function  or  cumulative  probabilities  are  computed  in 
terms  of  "error  functions,"  available  in  tabular  fashion  or  computer  routines. 
The  practical  application  of  this  formula  requires  only  that  we  substitute 
expressions  (5)  and  (6)  for  the  expectation  and  variance  in  (8)  to  compute  the 
distribution  of  any  resource  demand,  Dj.  The  cumulative  distributions  obtained 
in  this  way  are  generally  smoother  and  will  describe  the  distributions  well  in 
the  median  region.  In  the  distribution  wings,  however,  there  may  be  serious 
errors  in  this  approach  and  it  is  therefore  only  reconmended  for  quick 
estimations. 
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SECTION  FOUR 
DATA  ASSEMBLY 


4.0  Objective. 

The  objective  of  this  section  is  to  present  the  data  used  in 
demonstrating  the  Port  Systems  Requirements  Prediction  Methodology. 

4.1  Scope. 

This  section  describes  the  data  used  to  generate  the  Goal  Product, 
sources  of  the  data,  incorporation  into  the  computer  program  data  files,  and 
includes  comments  to  clarify  potentially  obscure  points. 

4.2  Discussion. 

Demonstration  of  the  Port  Systems  Requirements  Prediction  Methodology 
necessitates  the  development  of  data  files  for  ship  populations,  demands  for 
resources,  homeports,  operating  characteristics,  etc.  This  information  does 
not  exist  in  a  convenient  form  and  must  be  assembled  from  a  variety  of 
sources.  The  accuracy  of  its  selection  and  subsequent  development  bears 
directly  on  the  correctness  of  the  output  of  the  Metnodology. 

The  specific  information  required  by  the  computer  program,  PSP  as 
described  in  Appendia  A,  includes:  a  breakdown  by  class  of  ship  homeport 
assignments,  the  percentage  of  time  spent  in  port  at  cold  iron,  the  type  and 
amount  of  resources  demanded,  and  estimates  of  future  homeport  assignments 
for  1985  and  the  period  1995-2000.  Representative  data  from  existing  documents 
were  used  to  describe  these  characteristics  for  33  ship  classes  at  two  ports 
San  Diego,  CA,  and  Norfolk,  VA.  OPNAVINST  3111. 14U  was  used  to  establish 
homeportinq  assignments.  The  number  of  ships  in  each  class  with  homeports  of 
either  San  Diego  or  Norfolk  are  listed  in  data  files  DAI  and  DA2  in  Table  4-1  - 
Summary  of  Data. 
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Estimates  of  the  time  spent  in  port  at  cold  iron  were  derived 
using  Table  4-2  -  Ship  Operational  Modes,  obtained  from  reference  (b). 

The  percentage  of  time  spent  in  port  and  at  cold  iron  was  calculated  by 
averaging  the  percent  time  in  port  for  each  period  of  time,  the  percent 
time  at  cold  iron  for  each  period  of  time,  and  multiplying.  These  data  are 
listed  in  Table  4-1,  column  2. 

The  distribution  of  resources  was  restricted  to  electrical  demand. 

This  information  was  obtained  from  DM-25  using  KVA  values  for  Ships 
Service  (Hotel).  For  each  ship  class  listed,  a  single  value  of  KVA 
was  selected  to  represent  typical  electrical  demand  for  that  class  of  ships. 

The  electrical  demands  are  listed  in  column  three  of  Table  4-1  and  were  held 
constant  for  all  scenarios  under  study.  Estimates  of  the  1985  ship  class 
populations  were  based  on  Table  4-3,  extracted  from  reference  (d),  and 
involved  scaling  the  1981  ship  class  population  figures  up  or  down  accordingly. 
This  information  is  contained  in  data  files  DA3  and  DA4. 

The  ship  class  populations  in  data  files  DA5  and  DA6  are  estimates 
for  the  period  1995-2000  and  are  based  on  reference  (e).  Actual  figures  from 
reference  (e)  were  not  used  so  that  this  report  may  remain  unclassified. 

For  those  ship  classes  homeported  in  San  Diego  or  Norfolk  where 
information  on  resource  demands  and/or  percentage  of  time  in  port  at  cold 
iron  was  not  available,  an  estimate  was  made  based  on  ships  of  similar  type. 
Five  ship  classes  required  this  procedure  for  determining  their  percentage 
of  time  in  port  at  cold  iron.  Seven  classes  required  estimates  of  KVA 
demand. 
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SHIP 

TYPE 

PROBABILITY 

HOMEPC 

RT _ 

OF  DEMAND 

KVA 

19! 

31 

1985 

1995- 

■2000 

SD 

N 

N 

SD 

N 

AD 

.625 

3 

3 

3 

3 

2 

2 

AFS 

.381 

0 

3 

0 

3 

2 

5 

AO 

.500 

800 

0 

3 

4 

7 

0 

3 

AOE 

.314 

■  1  1 1  1 

0 

2 

0 

0 

0 

0 

AOR 

.349 

0 

2 

0 

2 

0 

2 

AR 

.713 

1714 

2 

1 

2 

1 

2 

1 

AS 

.352 

1250  SO 

2500  N 

2 

1 

3 

2 

2 

1 

ASR 

.569 

495 

2 

2 

1 

1 

2 

2 

ATF 

.657 

400  s 

1 

0 

1 

0 

1 

0 

CG 

.516 

2865 

9 

2 

11 

4 

12 

3 

CGN 

.388 

3293 

3 

5 

3 

5 

4 

6 

CV2 

.347 

6850 

3 

3 

3 

3 

6 

6 

CVN2 

.344 

11500 

0 

2 

0 

2 

0 

2 

CVT2 

.515 

935 

0 

1 

1 

0 

0 

DD 

.562 

1100 

9 

12 

13 

5 

5 

DDG 

.570 

1400 

9 

9 

13 

13 

5 

5 

FF 

.576 

1355 

12 

9 

15 

12 

24 

18 

FFG 

.599 

1012 

3 

1 

8 

6 

6 

3 

LCC 

.584 

311 

1 

1 

1 

1 

1 

1 

LKA 

.563 

1246 

3 

2 

3 

2 

3 

2 

LPA 

.566 

3000s 

0 

1 

0 

1 

0 

1 

LPD 

.516 

1636 

6 

7 

7 

8 

6 

7 

LPH 

.493 

4000 

3 

4 

3 

4 

3 

4 

LSD 

.576 

3000s 

7 

0 

7 

0 

7 

0 

LST 

.547 

3000s 

8 

0 

8 

0 

8 

0 

MSO 

.644 

3000  s 

2 

0 

2 

0 

2 

0 

SS14 

.701 

75 

3 

0 

3 

0 

3 

0 

SSN5 

.694 

1000 

12 

13 

22 

23 

12 

13 

AGFF 

.651 

467 

0 

1 

0 

1 

0 

1 

AGDS 

.651 

467* 

1 

0 

1 

0 

1 

0 

AGF 

.651 

467 

0 

1 

0 

0 

0 

1 

LHA 

.551 

3740 

2 

1 

3 

2 

2 

1 

PHM 

.651 

1006 

1 

0 

1 

0 

1 

0 

TOTAL 

33 

107 

89 

141 

120 

122 

95 

DATA  FILE: 

DAI 

DA2 

DA5 

0A6 

Probability  of  demand  estimated  due  to  lack  of  available  data. 

Aircraft  carriers  were  considered  as  two  groups  -  nuclear,  and 
non-nuclear.  The  training  carrier  CUT  was  considered  separately. 

electrical  demand  figures  are  for  ships'  service  (hotel)  only. 
They  do  not  Include  the  additional  power  required  for  Industrial 
work  load  operations  aboard  ship  when  In  for  repairs  or  overhaul. 

Includes  SS,  AGSS.  $ST,  and  LPS5  classes. 

SSN  and  SSBN  classes  grouped  together. 

Resource  demand  estimates  based  on  ships  of  similar  type. 


TABLE  4-1.  SUMMARY  OF  DATA 
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TABLE  4-2.  SHIP  OPERATIONAL  MODES 


TYPE 

FY-77 

INCREASE 

DECREASE 

AD 

3 

3 

_ 

AS 

1 

0 

- 

1 

AFS 

3 

3 

- 

- 

AGF 

1 

0 

- 

1 

AGFF 

1 

1 

- 

— 

AO 

4 

8 

4 

- 

AOE 

2 

0 

- 

2 

AOR 

2 

2 

- 

- 

AR 

1 

1 

- 

- 

AS 

1 

2 

1 

- 

ASR 

2 

1 

- 

1 

CG 

3 

5 

2 

- 

CGN 

4 

4 

- 

- 

CV 

3 

3 

- 

- 

CVN 

2 

2 

- 

- 

CW 

0 

1 

1 

- 

DD 

R 

11 

3 

- 

DDG 

9 

13 

4 

- 

FF 

q 

1? 

3 

- 

FFG 

1 

6 

5 

- 

LCC 

1 

1 

- 

- 

LHA 

1 

2 

1 

- 

LKA 

2 

2 

- 

- 

LPS  (NR) 

1 

1 

- 

- 

LPD 

6 

7 

1 

- 

LPH 

4 

4 

— 

- 

SSN 

12 

22 

10 

- 

TOTAL 

87 

117 

35 

5 

TABLE  4-3.  HOMEPORTED  SHIPS 


4.5 


SECTION  FIVE 


APPLICATION 

5.0  Objective. 

This  section  describes  the  application  of  the  Port  Systems  Requirement 
Prediction  Methodology  to  the  analysis  of  a  specific  utility  requirement  problem. 

It  emphasizes  the  procedures  necessary  to  acquire  a  goal  product  using  the  theory 
contained  in  Section  Three  and  the  data  assembled  in  Section  Four.  The  reader 
is  reminded  that  this  sample  problem  is  only  intended  to  demonstrate  the  Port 
Systems  Requirements  Prediction  Methodology.  Actual  application  of  the  Methodology 
requires  refinement  of  the  input  data. 

5.1  Scope. 

This  section  provides  detailed  information  concerning  the  computer 
procedures  that  are  necessary  to  implement  the  Methodology.  Typical  Methodology 
products  are  presented  with  a  discussion  of  the  information  that  can  be  derived 
from  them.  Finally,  recommendations  are  made  for  data  refinements  that  would 
be  desirable  for  further  application  of  the  Methodology. 

5.2  Data  and  Logic  Flow  Structure. 

A  schematic  diagram  of  the  data  flow  for  the  PSP  program  system  is 
presented  in  Figure  5-1.  The  input  data  set  is  prepared  on  disk  and  fed  to 
the  program  PSP  which  computes  the  means  and  variances  of  the  demands  and  the 
related  probability  distributions.  The  outputs  of  the  computations  are  distributed 
to  three  disk  files: 

1)  PSP. OUT,  which  is  an  image  of  the  hardcopy  output  of  the  program 
and  may  be  directed  either  to  a  line  printer  or  any  computer 
terminal . 
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PLOT  NO. 
1 
2 
3 


TYPE 

CUMULATIVE  PROBABILITY  VS  (LINEAR)  RESOURCE  AMOUNT 

CUMULATIVE  PROBABILITY  VS  (LOG)  RESOURCE  AMOUNT 

MEAN  AND  STANDARD  DEVIATION  OF  RESOURCE  USE  BY  SHIP  CLASS 


FIGURE  5-1  Schematic  Diagram  of  PSP  Code  Data  Flow 
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2)  PSP1.ARC,  which  is  an  archive  file  to  which  is  added  an  extract 
from  the  PSP. OUT  data  containing  just  the  probability  distributions 
necessary  to  produce  the  outputs  of  the  type  of  Figures  5-5 

and  5-6.  The  data  in  the  PSP1.ARC  file  is  reformatted  into 
plotter  compatible  format  using  the  program,  PLTFI1. 

3)  PSP2.ARC,  which  is  another  archive  file  to  which  we  add  an  extract 
of  the  PSP. OUT  data  containing  the  mean  and  variance  results. 

These  data  are  used  to  produce  a  plot  of  the  type  of  Figure  5-7. 

The  data  in  the  PSP2.ARC  file  are  reformatted  for  the  plotter 
using  the  program,  PLTFI2. 

Figure  5-2  presents  the  logic  flow  or  sequence  of  operations  for  the 
program  PSP1.  The  first  step  is  to  read  the  input  data.  This  may  be  accessed 
from  a  disk  file  when  constructed  off-line,  or,  in  another  version  of  the  program, 
it  may  be  entered  directly  at  the  keyboard  as  requested  in  conversational  mode 
with  the  computer.  The  first  calculation  is  to  obtain  the  probability  distribu¬ 
tion  of  equation  (1)  for  the  number  of  ships  present  in  each  class.  The  next 
computational  step  is  to  estimate  the  mean  numbers  of  ships  present  in  each  class 
as  given  in  equation  (3)  and  the  associated  variances  given  in  equation  (4). 

These  are  used  to  compute  the  means  and  variances  of  the  resource  demands  as 
specified  in  equations  (5)  and  (6). 

The  means  and  variances  of  the  resource  demands  are  fed  to  a  subroutine 
which  computes  the  distribution  of  total  demands  based  on  the  assumption  of  the 
validity  of  the  Central  Limit  Theorem,  i.e.,  that  the  total  demands  have  a 
Gaussian  distribution.  These  distributions  are  printed  at  this  time.  Finally, 
the  simulation  is  exercised  to  generate  the  desired  number  of  realizations  of 
the  numbers  of  ships  present  in  each  class  and  assess  the  total  resource  demands 
for  all  the  realizations  generated.  The  output  of  this  part  of  the  computation 
is  given  as  a  percentile  table  for  the  total  demand  for  each  resource. 
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READ  INPUTS 
FROM  TERMINAL 
OR  DISK 


FIGURE  5-2  Logic  Flow  Diagram  for  Program  PSP 
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If  more  inputs  are  indicated  either  through  terminal  interaction  or 
through  the  disk  data  set,  then  the  program  returns  to  the  data  reading  step. 
Otherwise,  execution  is  terminated  at  this  point. 

5.3  Data  Entry. 

Figure  5-3  presents  the  prompts  given  to  the  user  when  entering  a 
specific  example  of  a  data  set.  In  this  instance  the  program  was  reading  data 
from  a  disk  file  so  that  the  user  entries  do  not  appear  on  the  page. 


ENTER  NO.  OF  SHIP  CLASSES  AND  NO.  OF  RESOURCES  (213) 

ENTER  SCENARIO  DESCRIPTION  IN  72  CHARACTERS  OR  LESS 
ENTER  3  RESOURCE  NAMES  (NA4) 

ENTER  CLASS  NAME,  PROB.  OF  DEMAND  AND  RESOURCE  DEMANDS  (A4,  NF7.3) 
1 
2 


32 

33 

ENTER  SHIP  CLASS  POPULATIONS  (2413) 

ENTER  NO.  OF  SIMULATIONS,  NO.  DEPT,  AND  INITIAL  R.N.  (316) 
FIGURE  5-3  User  Prompts  for  Data  Entry 


An  initial  record  specifies  a  six  character  user  identi fication  code. 

The  next  line  of  data  gives  the  number  of  ship  classes  and  number  of  resources 
in  the  format  (213).  This  means  that  two  right  justified  three-digit  integers 
must  be  specified.  The  second  line  contains  up  to  72  characters  of  a  scenario 
description,  e.g.,  the  name  of  the  port  and  a  time  period  or  date.  The  names  of 
the  resources  must  then  be  specified,  each  having  up  to  four  characters 
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(alphabetic  and/or  numerical).  This  completes  the  header  record  set. 

The  main  body  of  the  input  data  is  the  resource  demand  matrix  giving 
the  demand  tor  each  resource  by  a  ship  in  each  class.  The  data  for  each  ship 
class  is  given  in  a  single  line.  The  first  four  characters  give  the  ship  class 
designation.  These  are  followed  by  the  expected  demands  for  the  resources  specified, 
each  demand  being  allowed  seven  characters  including  the  decimal  point. 

Two  closing  records  specify  the  ship  class  populations  and  the  simulation 
control  parameters.  The  number  of  ships  that  may  be  present  in  each  class  are 
given  in  the  first  record,  each  being  specified  as  a  right-justified  three-digit 
integer.  These  are  normally  the  numbers  of  ships  in  the  various  classes  that 
are  homeported  at  the  site.  The  last  record  specifies  the  number  of  simulations 
to  be  performed  (zero  indicates  no  simulation  desired),  the  number  of  (largest) 
demands  to  be  retained  for  developing  percentile  tables  for  each  resource,  and 
the  initial  random  number  which  can  be  any  odd  integer.  Each  of  the  numbers  in 
the  last  record  is  given  as  a  right-justified  six-digit  integer. 

The  computer  returns  a  page  presenting  a  formatted  description  of  the 
input  as  illustrated  in  Figure  5-4.  This  may  be  reviewed  by  the  user  for  correct¬ 
ing  errors  and  filed  for  record  keeping  purposes. 

5.4  Practical  System  Operation. 

This  PSP  model  is  best  executed  in  time-sharing  mode  using  a  computer 
terminal.  It  was  developed  on  the  ADP  network  and  runs  on  a  DEC  10  computer. 

The  graphics  output  routines  are  written  for  a  Tectronix  4051  computer/plotter. 

While  there  are  some  elements  of  the  PSP  program  that  are  system  specific,  it  is 
readily  modified  to  work  on  any  computer.  The  graphics  routines  will  always 
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NORFOLK,  1995-2000 

NUMBER  OF  CLASSES*  33  NUMBER  OF  RESOURCES*  3 

NUMBER  OF  REALIZATIONS  REQUESTED=  100 
NUMBER  KEPT*  100 


INITIAL  RANDOM  NUMBER*  *190373 
RESOURCE  DEMAND  MATRIX  IS  AS  FQL.LOUS : 


KUA 

KWDC 

NQ  . 

AD 

2200 . 000 

0.000 

1.000 

AFS 

1184.000 

0.000 

1.000 

AO 

800 . 000 

0.000 

1.000 

AOE 

1090 .000 

0  .  000 

1.000 

AOR 

1012.000 

0.000 

1.000 

AR 

1714.000 

0.000 

1  .  000 

AS 

2500.000 

0.000 

1.00  0 

ASR 

495.000 

0.000 

1.000 

ATF 

400.000 

0.000 

1.000 

CG 

2865.000 

0.000 

1,000 

CGN 

3293.000 

0  .  000 

1.000 

CV 

6850 .000 

0  ,  000 

1  .000 

CUN 

11500.000 

0.000 

1.000 

CUT 

935.000 

0.000 

1.000 

DD 

1100.000 

0.000 

1.000 

DDG 

1400.000 

0  .  000 

1.000 

FF 

1355.000 

0.000 

1.000 

FFG 

1012.000 

0  .  000 

1  .  000 

(..CC 

31  1  .000 

576.000 

1  .  000 

I..K  A 

1246.000 

0 .000 

1.000 

LPA 

3000.000 

0 .000 

1.000 

LPD 

1636.000 

0 .000 

1.000 

L.PH 

4000.000 

0  .000 

1.000 

LSD 

3000.000 

0.000 

1  .  000 

L  ST 

3000.000 

0.000 

1.000 

MSO 

3000.000 

0.000 

1.000 

SS 

75.000 

165.000 

1.000 

SSN 

1000.000 

0.000 

1.000 

AGFF 

467.000  . 

0.000 

1.000 

AGDS 

467.000 

0  .000 

1.000 

AGF 

467.000 

0.000 

1.000 

LHA 

3740 . 000 

0.000 

1.000 

PHM 

100. 000 

0  .000 

1.000 

PROB 

0 . 625 

n 

«_ 

0 .381 

5 

0 . 500 

3 

0.314 

0 

0 . 349 

') 

L.. 

0.713 

1 

0.352 

1 

0.569 

«~ 

0 . 657 

0 

0.516 

3 

0 . 388 

6 

0.347 

6 

0 . 344 

2 

0.515 

0 

0 . 562 

5 

0.570 

5 

0 .576 

18 

0.599 

3 

0.594 

1 

0.563 

2 

0.566 

1 

0.516 

7 

0 . 493 

4 

0.576 

0 

0.547 

0 

0.644 

0 

0.701 

0 

0.694 

13 

0.650 

1 

0.650 

0 

0 . 650 

1 

0.550 

1 

0 .650 

0 

FIGURE  5-4  Sample  Formatted  Page  Returned  for  Operator  Review 
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require  a  Tektronix  system,  however,  as  new  programs  would  have  to  be  generated 
to  work  on  another  system. 

When  operating  with  a  real  time  input  version,  the  user  merely  responds 
to  terminal  requests  for  input  once  the  six  character  run  identifier  has  been 
given.  For  construction  of  the  disk  data  file,  any  system  text  editor  can 
be  used. 


5.5  Methodology  Output. 

Section  2  introduced  the  goal  product  of  the  Port  Systems  Requirements 
Prediction  Methodology,  i.e.,  a  curve  depicting  the  probability  of  exceeding 
the  demand  for  a  selected  utility  for  a  given  port  condition  and  discussed  its 
utilization.  The  computer  programs  contained  in  Appendix  A  provide  this  product 
as  an  automatically  graphed  output.  Examples  of  these  goal  products  are  given 
in  Figures  5-5,  5-6,  and  5-7.  These  curves  present  unique  ways  in  which  the  primary 
output  of  the  Port  System  Requirements  Prediction  Methodology  may  be  used. 

These  techniques  are  described  as  follows: 

Figure  5-5.  Three  sets  of  data  are  represented  here.  They  are  the 
ships  homeported  in  San  Diego  in  1981,  the  ships  that  would  be  in  San  Diego  in 
1985  if  the  fleet  were  to  expand  in  size  comparable  to  that  depicted  in 
Sewell's  Point  Master  Plan  (1985),  and  a  1995  -  2000  fleet  determined  by  as¬ 
suming  that  the  projections  of  reference  (e)  are  realistic.  Applying  the 
technique  of  Section  2,  it  can  be  seen  that  a  twenty  percent  probability  of 
the  demand  for  electricity  exceeding  the  supply  would  grow  to  70%  and  90%  for 
the  1995  and  1985  Scenarios,  respectively,  if  there  is  no  increase  in  supply. 
Alternatively,  to  maintain  the  probability  of  failing  to  meet  demand,  the 
supply  must  be  increased  16-24000  KVA.  This  type  of  analysis  can  be  used  to 
determine  the  impact  of  changes  to  current  homeport  assignments,  base  closures, 
or  to  determine  the  adequacy  of  utility  supply  for  alternative  concepts  of  future 
fleets.  Other  applications  exist  at  the  pier  level.  For  example,  utility 


adequacy  for  a  specific  pier  loading  could  be  established  and  compared  to 
alternative  pier  loading  plans. 

Figure  5-6.  The  output  depicted  in  Figure  5-6  is  formatted  the  same 
as  Figure  5-5.  However,  the  manner  in  which  the  initial  conditions  are  established 
makes  this  a  unique  output.  Here  the  per  ship  current  demand  is  incremented  by 
a  predetermined  amount  and  the  resultant  increase  in  failure  probability  can 
be  observed.  Using  this  technique,  a  port  planner  could  determine  when 
improvements  to  ships,  installed  in  accordance  with  the  Fleet  Modernization 
Program,  would  aggregate  to  the  point  of  affecting  a  particular  pier.  However, 
certain  data  files  would  have  to  be  constructed  before  this  could  be  achieved. 

This  would  be  expensive  on  a  Navy-wide  basis,  but  it  could  be  very  effective 
at  selected  ports  where  the  use  of  selected  facilities  could  be  statistically 
characterized. 


Both  Figures  5-5  and  5-6  result  from  applying  the  Central  Limit 
Theorem  (solid  line)  and  the  Simulator  (points)  for  the  different  scenarios 
under  study. 

Figure  5-7.  Figure  5-7  presents  a  bar  chart  indicating  the  mean  and 
standard  deviation  of  the  usage  of  the  given  resource  determined  by  the  popula¬ 
tion  and  deployment  cycle  of  the  different  ship  classes.  This  chart  enables  one 
to  identify  the  critical  classes  for  consumption  of  the  given  resource.  In 
this  regard,  and  referring  to  the  use  of  Figure  5-6,  this  chart  could  provide 
the  initial  guidance  for  determining  the  key  classes  of  ships  that  one  should 
include  in  a  review  of  the  Fleet  Modernization  Program  as  a  prerequisite  for 
studying  selected  piers. 


5.6  Conclusions  and  Recommendations. 

All  that  has  been  written  in  this  report  up  to  here  has  been  directed  at 
establishing  an  algorithm,  or  methodology,  by  which  statistical  techniques  might 
be  used  for  port  analysis.  A  set  of  data  has  been  operated  upon  by  the  programs 
developed  for  the  PSRPM  and  the  automatic  generation  of  output  has  been  demonstrated. 
The  shortcomings  of  the  process  are  those  that  were  built  in  by  the  initio!  data 
used.  However,  the  important  fact  to  note  is  that  the  methodology  worked  and 
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FIGURE  5-6  PROBABILITY  OF  EXCEEDING  SPECIFIED 
DEMAND  FOR  KVA. 
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FIGURE  5-7  MEAN  and  STANDARD  DEVIATION  of  demand  for  Electrical  Power. 


since  it  did,  a  significantly  powerful  and  documented  tool  is  now  available  to 
port  planners  and  designers.  The  challenge  is  to  learn  to  use  the  tool  and 
to  describe  problems  that  can  use  its  approach.  In  this  sense,  some  initial 
recommendations  concerning  data  bases  to  establish  and  possible  near  term  studies 
to  pursue  for  validating  or  refining  the  Port  Systems  Requirements  Prediction 
Methodology  can  be  made.  These  are: 

a.  Data  Bases 

(1)  An  immediate  refinement  to  the  Methodology  can  be  made  by 
using  ship  class  electrical  demand  data  that  is  based  on  actual 
measurements. 

(2)  Deployment  cycle  and  turnaround/workup  schedule  information 
can  be  used  to  establish  the  amount  of  time  a  ship  is  operating 
out  of  its  home  port.  This  will  be  used  to  establish  the  total 
demand  for  a  utility  that  will  be  made  by  selected  ship  groupings. 

(3)  Various  conceptions  of  total  fleet  composition  can  be 
developed  based  on  alternative  strategies  of  fleet  development. 

For  example, 

A  Type  I  Fleet  could  be  a  projection  of  current  types  and 
levels  to  some  maximurg, 

A  Type  II  Fleet  could  be  characterized  by  a  concentration 
of  capital  ships,  i.e.,  CVAs,  CGs,  Aegis  ship,  etc., 

A  Type  III  Fleet  could  be  characterized  by  high  concentration 
of  small  combatants,  and 

A  Type  IV  Fleet  could  be  the  best  estimate  of  a  1990  Navy. 

These  fleets  could  then  be  used  in  selected  problems 
and  would  serve  as  a  basis  for  comparing  port  concepts. 

b.  Near  Term  Studies 

An  immediate  objective  of  near  term  work  with  the  Port  Method¬ 
ology  is  its  validation  and/or  refinement.  This  could  be  achieved  by 
using  electrical  measurement  data  to  establish  per  ship  demand  and 
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then  using  this  data  file  to  predict  demand  relationships  at 
the  pier  level.  The  output  of  such  a  study  could  be  used  to 
assess  the  impact  of  selected  MILCON  Projects  and/or  to  serve 
as  a  basis  of  an  acquisition  strategy.  Finally,  as  familiar¬ 
ity  with  the  process  is  gained,  procedures  could  be  developed 
for  presenting  NAVFAC  and  Fleet  planners  with  an  analytical 
review  of  existing  and  planned  facilities. 

Additional  studies  certainly  can  be  developed,  but  at 
this  stage  of  the  development  of  the  Port  Systems  Requirements 
Prediction  Methodology,  it  would  be  premature  to  describe  them 
Rather,  the  development  and  execution  of  a  meaningful  initial 
study  using  measured  data  should  be  the  primary  objective.  In 
this  regard,  any  on-going  or  near-term  utility  measurement 
project  could  be  used  for  this  purpose. 
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APPENDIX  A 


PSP  COMPUTER  PROCESSING  DOCUMENTATION 

PROGRAMS 

PSP  -  PSP. FOR  is  the  main  processing  program  which  computes  the  central  limit 
theorem  estimate  and  the  Monte  Carlo  Simulations  of  port  system  demand  for 
specified  port  resources  and  predicting  scenarios.  On  each  run  of  PSP. FOR 
(i.e.  each  specified  scenario)  two  sets  of  data  are  written  to  an  archive  storage 
file.  The  first  set  of  data  output  is  the  port  system  demand  profile  which  con¬ 
tains  both  the  central  limit  theorem  and  the  simulation  projections.  This 
archive  file  is  named  PSP! .ARC  and  resides  on  the  system  hard  disc.  The  second 
set  of  data  output  is  the  ship  class  demand  profile  which  contains  the  computed 
mean  and  standard  deviation  of  resource  demand  for  each  ship  class.  This 
archive  file  is  named  PSP2.ARC  and  also  resides  on  the  system  hard  disc.  Both 
PSP.l.ARC  and  PSP2.ARC  are  sequential  files  using  operator  defined  five  letter 
scenario  codes  to  identify  each  scenario  input  to  PSP. FOR. 


A-l 


00010 
>00020 
00030 
00040 
00050 
00060 
00070 
OOOSO 
00090 
00100 
00110 
00120 
00130 
00140 
00150 
00160 
00170 
00180 
00190 
00200 
00210 
00220 
00230 
00240 
00250 
00260 
00270 
, 00280 
00290 
00300 
00310 
„  00320 
00330 
•  00340 
’  00350 
00360 
00370 
00380 
00390 
00400 
00410 
00420 
00430 
00440 
00450 
00460 
00470 
00480 
00490 
00500 
00510 


C 

C 

C 

C 


C 

C 

C 

C 


c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 


THIS  PROGRAM  SIMULATES  PORI  SERVICE  DEMAND  OVER  A  NUMBER  Or  TIME 
PERIODS  AND  DETERMINES  THE  PROBABILITY  DISTRIBUTION  OF  THE  DEMAND 
FOR  A  NUMBER  OF  RESOURCES. 

READ  THE  INPUT 

OPEN (UNIT *5 > ACCESS^ ' SEQIN ' ,FILE='PSP . DAT ' , DEVICE* ' DSK ' > 
0PEN(UNIT=8  > ACCESS- ' APPEND ' ,FILE~'PSP2 . ARC' , DEVICE- ' DSK ' > 

OPEN  ( UNIT  =  6 ,  ACCESS^ '  SEQOUT  '  ,FILE*'PSP  .OUT'  ,  DEVICE-- '  DSK  '  ) 

1  CALL  DESCR 

COMPUTE  MEANS  AND  VARIANCES  OF  CLASS  DEMANDS 
CALL  STATS 

COMPUTE  THE  CUMULATIVE  PROBABILITY  DISTRIBUTIONS 
CALL  INIT 

PERFORM  SIMULATION 
CALL  SIMUL 
OUTPUT  THE  RESULTS 
CALL  OUTPUT 

CLOSE  (UN  IT  =7 , ACCESS- ' APPEND '  ,FILF.~'PSP  .ARC'  ,  DEVICE* '  DSK. '  ) 

CLOSE  (UNIT55  6  ,  ACCESS-  '  SEQOUT  '  ,FILE*'PSP  .OUT'  ,  DEVICE* '  DSK  '  ) 

CLOSE (UNIT «8 ,  ACCESS"  '  APPEND  '  ,FILE*'PSP2. ARC'  ,  DEVICE* '  DSK  '  > 


END 

SUBROUTINE  DESCR 


THIS  ROUTINE  INPUTS  THE  REQUIRED  DATA 
THE  INPUT  PARAMETERS  APE  AS  FOLLOW i 

NCLASS  =  NUMBER  OF  SHIP  CLASSES  SERVICED  BY  FACILITY 
=  0  TO  STOP 

=  NEGATIVE  If  ONLY  THE  SHIP  POPULATION  CHANGES  FROM 
PREVIOUS  CASE  RUN 


NUMBER  OF  RESOURCES  OF  INTEREST 
1  FOR  A  DUMP  OF  THE  ’  ">UT  <0  OTHERWISE) 

NAME  OF  ITH  RESOUCE 
NAME  OF  ITH  SHIP  CLASS 

PROBABILITY  THAT  MEMBER  OF  ITH  CLASS  REQUIRES  SERVICE 
AMOUNT  OF  JTH  RESOURCE  REQUIRED  BY  ITH  SHIP  CLASS 
TOTAL  RELEVANT  POPULATION  IN  ITH  SHIP  CLASS 
NUMBER  OF  REALIZATIONS  TO  BE  GENERATED 
NUMBER  OF  HIGH  DEMANDS  OF  EACH  RESOUCE  TO  BE  RETAINED 
INITIAL  RANDOM  NUMBER  (=0  FOR  STANDARD  R.N.  SET) 

MCQDE  *  6  CHARACTER  ID.  CODE 
DIMENSION  MCODE(fe) 

CQMMON/DATA/NCLASS , NRS ,NPER , NKE , IX , NAMER ( 20  > , NSHIP  <  50  > , PR  <  50  > , 
1NAMEC (50  > ,R ( 20 ,50 ) , I DUMP ,MSG(1B) 

OPEN  <  UNIT«7 , ACCESS* ' APPEND ' , F I LE  ~ 'PSP . ARC ' , DEVICE* ' DSK ' > 


NRS 

IDUMP  * 
NAMER  ( I ) 35 
NAMEC ( I ) » 
PR  <  I ) 

R  ( J  ,  I )  = 

NSHIP ( I > - 

NPER 

NKE 

IXX 


IX=4973127 

NCL*NCLASS 

READ <5, 105,END*5>  MCODE 


105  FORMAT (6A1 ) 

WRITE (7, 105)  MCODE 

WR ITE (8 , 1 05 )  MCODE 

READ (5,100) NCLASS ,  NRS ,  I DUMP 

IF (NCLASS . EQ . 0 )GQ  TO  5 

READ(5,89)MSG 

A- 2 


J)  052  0 
00530 

0  0540 
■10550 
*0  0560 
00570 

106 

110580 
130590 
00600 
iQOol  0 

999 

'00620 

00630 

2 

00640 

10650 

3 

0  0660 

4 

00670 

00680 

10690 

0  0700 

1  07 

0  071  0 
00720 

0  0730 

5 

00740 

00750 

00760 

88 

0  0  7'7  0 

89 

0  0780 

90 

00790 

91 

00800 

92 

00810 

00820 

93 

00830 

94 

0  0840 

95 

0  0350 

1  0  0 

0  0860 

101 

0  0  870 

1  0  2 

0  0880 

1  0  3 

0  0890 
00900 
0  0  91  0 
0  0920 
0  0930 
0  0  940 
0  0930 
0  0960 
0  0970 
00980 


UR ITE (7 , 89 )  HSG 
WRITE <8, 89)  HSG 
IF (NCLASS . LT . 0 )GO  TO  3 
UR ITE( 7 ,106)  NCLASS , NRS 
WR  ITE  <  8 , 1  06 )  NCLASS , NRS 
F0RMAT(I3,3X,I3) 

READ (5,101) (NAMER (I) ,1=1  ,  NRS) 

WR ITE ( 7 , 999 )  <  NAHER ( I ) ,1=1 ,NRS) 

UR ITE (8, 999)  ( NAHER  < I ) ,1=1 ,NRS) 

FOR HAT (3<A4,3X) ) 

DO  2  1=1 , NCLASS 

RE AD <5, 102) NAHEC ( I ) ,PR ( I ) , <R < J , I ) , J= 1 ,NRS) 

GO  TO  4 

NCLASS=NCL 

CONTINUE 

READ<5, 100) (NSHIP (I ) ,1=1 , NCLASS) 

READ  (5, 1Q3)NP£R  ,  NKE  ,  I. XX 
WRITE <7,107)  NPER , NKE 
FOR HAT  < 16 , 3X , 16 ) 

IF C I XX .NE . 0 ) IX=IXX 
GO  TO  6 

CLOSE  <  UNIT=5 , ACCESS= ' SEQ1N ' ,FILE='PSP .DAT' ,DEVICE='DSK ' ) 

CLOSE ( UNIT=7 , ACCESS = 'APPEND' ,FILE='PSP .ARC' ,DEVICE='DSK ' ) 

STOP 

FORMAT  <  '  ENTER  SENARIO  DESCRIPTION  IN  72  CHARACTERS  OR  LESS') 
FORMAT < 1 3A4 ) 

FORMAT (///'0ENTER  NO.  OF  SHIP  CLASSES  AND  NO.  OF  RESOURCES < 213 ) 


'  ) 

OURCE  DEMANDS (A4 


SHIP  CLASS  P0PULATI0NS(24I3) ' ) 
NO.  OF  SIMULATIONS, NO.  KEPT, AND 


3) 


INITIAL  R . N .  (316) '  ) 


C 

C 


FORMAT (  '  ENTER', 14,'  RESOURCE  NAMES  ( NA4 ) ' ) 

FORMAT ( '  ENTER  CLASS  NAME,  PROP.  OF  DEMAND  AND 
1 , NF7 . 3 )  ’  ) 

FORMAT ( '  ',13) 

FORMAT ('  ENTER 
FORMAT  < '  ENTER 
FORMAT (2413) 

FORMAT (20A4) 

FOR MAT <A4 , 1QF7 
FORMAT (316) 

,  RETURN 
END 

SUBROUTINE  I NIT 

THIS  ROUTINE  INITIALIZES  THE  CUMULATIVE  DISTRIBUTION  FUNCTION 
COMPUTED  ON  THE  BASIS  OF  THE  BINOMIAL  DISTRIBUTION. 

COMMON 'COMP /PR OB ( 1000) ,NST(51 ) 

CQMMON/AMER IT/BIG (20 , 1 00 ) , IPOS ( 20 , 1 0 0 ) , NKEPT , REQ ( 20 ) 
COMMON/DATA/NCLASS, NRS, NPER , NKE , I X , NAMER ( 20 ) ,NSHIP(50) ,PR (SO) , 
l NAHEC <50 ) ,R(20 ,50 ) , I DUMP , MSG ( 1 8 ) 

NST (  1  )= 0 


00990 
01000 
01010 
01020 
01030 
01040 
01050 
01  060 
01070 
oi  oeo 

01090 

01100 

OHIO 


DO  1  1=1,  NCLASS 
NN=NSHIP ( I ) 

FN=NN 
NN1 =NN+1 

NST ( I + 1 ) =NST ( I ) +NN1 
P=PR ( I  ) 

Q=1 . -P 

PQ--P  /Q 

K=NST(I)+1 

PP=Q**FN 

PROB ( K ) =PP 

IF ( NN . EQ . 0 ) GO  TO  1 

DO  1  J=1 , NN 


.1120 

01130 

11140 

.1150 

01160 

011.70 

n  iso 

J  1  I  9  0 
01200 
1121  0 
11220 
01230 
01240 
11250 
J  1  2  6  0 
01270 
1)  2 SO 
J  1  2  9  0 
01300 
11310 
i  1320 
0  1330 
01340 
11350 
ill  360 
01370 
11380 
.11390 
0  1400 
3  1  4 1  0 
11420 
01430 
0  I  44  0 
)  1  430 
I)  1  4  6 1 

0  i  471) 
1  i  49 0 
'.31490 
01500 
Or  1  ' :  3  0 
0  1520 
0  1530 
01540 
01530 
01560 
01570 
01580 
01590 
01  600 
0161  0 
01620 
0  1630 
01640 
01650 
0  1660 
01670 
0  1  68  0 
0  I  69  0 
0  1700 
0171  0 


K  =K  + 1 

PP=PP*PQxFLOAT <  NN1 - J ) /FLOAT  (  J ) 

PR0B(K)=PR0B(K-1 >+PP 

1  CONTINUE 
NKEPT=0 
RETURN 
END 

SUBROUTINE  SIMUL 

C  THIS  ROUTINE.  SIMULATES  PORT  SERVICE.  DEMAND  FOR  A  SET  OF  SPECIFIED 

C  RESOUCES. 

COMMON/COMP /PROBi 1000) ,NST<51 ) 

COMMON/ AMER IT/BIG ( 20,1 0  0  ) , IPOS (20 , 1 0  0 ) ,NKEPT , REQ (20 ) 

COMMON  /  DAT  A/NCL.  ASS  , NR S , NPER , NKE , I X ,  NAMER  ( 20  )  ,NSHIP(50>  ,PR(5G)  , 

1 NAMEC ( 50  ) ,R(20,50) ,  IDUMP ,MSG< 18) 

WRITE (6 , 1 00 )  (NAMER (II ) ,11=1, NR S) 

DO  1  1=1 , NPER 
DO  2  L-l.NRS 

2  REQ ( L ) =0 , 

DO  3  J= 1 , NCLASS 
CALL  R AN DU (IX,IY,XX) 

I X  =:  I Y 
NN=NST ( J ) 

NM=NSHIP ( J) 

DO  4  K--1  ,NM 

IF (XX ■ LT . PROBCK+NN) )GO  TO  5 

4  CONTINUE 

5  FN=FLOAT(K-l ) 

DO  6  L.  =  l  ,  NRS 

6  REQ ( L ) =REQ ( L ) +FN*R ( L , J  > 

3  CONTINUE 

UP ITE ( 6 , 1 0 1 )  I, (REQvL) ,L=1 ,NRS) 

ILOC-I 

1  CALL  PLACE(ILOC) 

100  FORMAT ( /  / '  SIMULATED  REALIZATIONS  OF  RESOURCE  DEMANDS'/ 

1'  ' ,5X,6(3X,A4,2X) >  I 

101  FORMAT (  '  ',14, 8F9 . 1  )  I 

RETURN  J 

END  *  J 

SUBROUTINE  OUTPUT  j 

C  THIS  ROUTINE  OUTPUTS  THE  LARGEST  DEMANDS  AND  THEIR  GENERATION  NO! 

COMMaN/COMP/PROBdOOD)  ,NST(51  )  J 

COMMON/AMER IT/B IG ( 20,1  00  ) , IPOS (20 , 1  00  )  NKEPT , REQ < 20 )  I 

C.OMMON/DATA/NCLASS  ,  NRS  ,  NPER  ,  NKE  ,  IX  ,  NAMC.s  (  20  )  ,  NSHIP  (50  >  ,PR  (50  )  ,  J 
1 NAMEC ( 50 ) , R ( 2  0 ,50) ,IDUMP ,MSG(18)  I 

WRITE (6, 100)  NKE , (NAMER ( I ), 1=1 , NRS)  I 

DO  2  J=t,NRS  I 

WRITE (7,102)  <  BIG ( J , I )  ,  1  =  1, NKE)  J 

2  CONTINUE  I 

DO  3  1*1, NKE  I 

WR ITE (6,101)  I  ,  (BIG( J, I ) , IPOS( J , I ) , J  =  1 ,NRS)  j 

3  CONTINUE  I 

102  FORMAT ( 20 ( 2X , El 1 . 4 ) )  j 

100  FORMAT ( ' THE  LARGEST ' , 14  ,  '  DEMANDS  FOR  EACH  RESOURCE  WERE  AS  FO  J 

1 LLOW  '  / '  ' , 6 ( 1  OX , A4 ) )  I 

101  FORMAT  ('  ' ,I3,5(F9, 1 ,14) )  I 

RETURN  I 

END  I 

SUBROUTINE  PLACE(ILOC)  I 

C  THIS  ROUTINE  CHECKS  TO  SEE.  IF  ANY  OF  THE  COMPUTED  RESOURCE  DE.M| 

C  IN  THE.  IL.OCTH  REALIZATION  LIE  IN  THE  HIGHEST  GROUP  I 


fl  720 
1730 

11740 
1750 
1760 
01770 

11780 
1790 
01800 
f  1810 
1.1820 
1)1830 
01840 
f  1850 
1)1360 
01870 
1US80 
1890 
019  0  0 
-01910 
1)1920 
0  1930 
01940 
'.195  0 
'  J  1  9  6  0 
01970 
D1980 
11990 
02000 
-02010 
12020 
02030 
02040 
12050 
J7060 
02070 
- 12080 
)  2090 
02100 
02110 
12120 
02130 
02140 
'  12150 
12160 
02170 
12180 
02190 
02200 
. 02210 
12220 
02230 
02240 
' 12250 
02260 
0227  0 
02280 
02290 
02300 
02310 


COMMON/COMP /PR OB ( 1000) ,NST (51 ) 

COMMON/ AMER IT/BIG ( 20 , 1 0  0  > , IP OS <20,100) ,  NKEPT , REQ ( 20 ) 
COMMON/DATA/NCLASS , NRS , NPER ,  NKE , IX , NAMER ( 20 ) ,NSHIP<50> ,PR<50> 
1NAMECC50) , R ( 20 , 50 ) , IDUMP , MSG (18) 

IF(NKEPT  .  GT  .  NKE)NKEPT=NKL" 

NK=NKEPT 
NK2-NK+2 
NK3  =  NK  +3 
NKEPT=NKEPT+1 
DO  12  1  =  1, NRS 
J=  1 

IF (NKEPT . EQ . 1 )GO  TO  2 

IF<  NKEPT . GT . NKE . AND . REQ  < I > . LT . BIG  < I , NK ) >GO  TO  12 
DO  3  J=1 , NK 

IF <REQ< I > . LE ,BIG(I,NKEPT-J) )GO  TO  4 

3  CONTINUE 

4  IF( J .EQ. 1 )GO  TO  2 

MOVE  UP  (J-l)  VALUES  AND  POSITION  INDICATORS 
DO  5  K  -  2 ,  J 

IF ( NK3-K . GT . NKE ) GO  TO  5 
BIG(I,NK3-K)=BIG(I,NK2~K) 

IPOS ( I , NK3-K ) =IPOS ( I ,NK2-K > 

5  CONTINUE 

INSERT  NEW  VALUES 
2  BIG < I , NK2- J >  =REQ ( I ) 

IP0S(I,NK2-J)=IL0C 
12  CONTINUE 
RETURN 
END 

SUBROUTINE  DUMP 

THIS  ROUTINE  DUMPS  THE  INPUT  DATA 

COMMON/DATA/NCLASS , NRS , NPER , NKE , IX , NAMER ( 20 ) ,NSHIP (50 ) ,PR (50 ) 
1 NAMEC ( 50 ) , R ( 20 , 50 ) , IDUMP ,MSG( 1 8) 

DATA  LAB/ ' PROB ' / 

WRITE (6, 99)  MSG 

WRITE (6,1 00 )NCLASS, NRS, NPER , NKE, IX, (NAMER ( I > , I = 1 , NR S ) ; LAE 
DC)  1  I  =  1  ,  NCLASS 

WRITE(6, 101 )  NAMPC ( I ) , ( R ( J , I ) , J= 1 , NR S ) , PR ( I ) 

1  WR ITE ( 6 , 1 02 )  NSHIP(I) 

99  FORMAT ( 1 8A4 ) 

100  FORMAT ( 'NUMBER  OF  CLASSES^ ' , 1 4 , 6X , ' NUMBER  OF  RESOURCES=  '  , 
114/' NUMBER  OF  REALIZATIONS  REQUESTED= ', 15/ ' NUMBER  KEPT*' 
2,14/'  INITIAL  RANDOM  NUMBERS  ,110/ 

3'  RESOURCE  DEMAND  MATRIX  IS  AS  FOLLOWS:'/'  ' , 3X , 6 ( 6X , A4 ) ) 

101  FORMAT (A4, 5F 10,3) 

102  FORMAT (65X, 15) 

RETURN 

END 

SUBROUTINE  STATS 

THIS  SUBROUTINE  COMPUTES  MEANS  AND  VARIANCES  OF  RESOURCE 
REQUIREMENTS  OF  THE  CLASSES  AND  COMPUTES  THE  PERCENTILES 
OF  THE  TOTAL  RESOURCE  USE  BY  THE  CENTRAL  LIMIT  THEOREM. 
COMMON/DATA/NCLASS, NRS, NPER , NKE , IX , NAMER (20 ) ,NSHIP (50 ) ,PR  <50  > 
1 NAMEC (50) ,R( 20,50) , IDUMP , MSG ( 1 8 ) 

DIMENSION  OUT <20, 2) , SUM (20) , VAR (20) , PERC ( 23 ) , PMULT ( 23 ) 

DATA  PERC/99 .  , 98 .  , 97 .  ,96 .  , 95 .  , 90 .  , 85 .  , 80 .  , 75 .  , 70 .  , 60 .  , 50 .  , 
140.  , 30.  ,  25.  , 20. , 15. ,10. ,5.  ,4,  ,3., 2.  ,1./ 

DATA  PMULT/-2 .327,-2 . 054,-1 .881,-1.751,-1 .645,-1 .382,-1 .037, 
1-0 . 842,-0 . 674,-0 . 524,-0 . 253,0 ,0,0. 253,0 .524 , 0 . 674,0 . 842,1 . 037 
21 . 382,1 .645,1 .751,1.881,2.054,2. 327/ 


|l2320 

DATA  NPERC/23/ 

it  2330 

WRITE (6,99)  MSG 

-02340 

WRITE <6, 100)  <  NAMER ( J ) , J=1 ,NR3> 

■02350 

DO  1  J--1  ,  NRS 

^02360 

SUM ( J ) =0  . 

02370 

1 

VAR < J)=0 . 

|02380 

DO  3  1=1 , NCLA5S 

•  02390 

AN=FLOAT (NSHIP ( I ) Jx-PRCI) 

0  240  0 

FN=AN*(i  .-PR  (I)  ) 

j0241  0 

DO  2  J=1,NRS 

(02420 

OUT  <  J , 1 )=R(J,I)XAN 

02430 

V=R ( J , I )  **2*FN 

02440 

OUT ( J  ,  2 ) =SQRT  <  V ) 

0  2450 

SUM < J ) =SUM ( J  >  +OUT ( J ,  1 ) 

02460 

2 

VAR  ( J  )  =VAR  ( J )  +V 

02470 

WRITE <6,101)  NAMEC ( I > , <  OUT ( J ,  1 ) , J  =  1 ,NR 

0  246  0 

WR ITE ( 6  j 102)  <OUT< J,2) , J=1 ,NRS> 

02490 

WRITE (8, 109)  <  OUT  <  J  ,  1 ) ,OUTvJ,2) ,J  =  1 ,NR 

02500 

109 

FORMAT ( 6 ( 2X  ,  F 1  0 . 3) ) 

02510 

3 

CONTINUE 

02520 

DO  4  J=1 ,NRS 

02530 

4 

VAR  (  J ) =SQRT ( VAR  <  J ) ) 

02540 

99 

FORMAT  < // '  ' , 1 8A4/ '  MEANS  AND  STANDARD 

0255  0 

1NDS ' ) 

02560 

100 

FORMAT  (  '  '  >4X,6U0X,A4)  ) 

0  257  0 

101 

FORMAT ( '  '  ,A4,6(4X,F10 , 3) ) 

0  2580 

1  02 

FORMAT  ('  '  ,2X,6<'  '^lO.S)) 

02590 

1  0  3 

FORMAT ( ' 0 ALL ' , 6 < 4X , F 1 0 . 3 ) ) 

02600 

105 

FORMAT < '  '  ,F6 . 1 ,4X,6(F10 .3,4X) ) 

02610 

DO  5  1=1 , NPERC 

02620 

DO  6  J=1 ,NRS 

02630 

6 

OUT  v  J  ,  1  )  =SUM  ( J )  +PMUI..T  ( I  >  *VAR  (  J  ) 

02640 

WRITE (7 ,96)  ( OUT <  J , 1 > ,J=1 ,NRS) 

02650 

5 

CONTINUE 

0  2  6  6  0 

96 

FORMAT  ( 20  (2X,EU  .4)) 

0  2670 

RETURN 

02680 

END 

DEMA* 
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PLTFI1 


-  PLTFI1.F0R  is  the  program  which  generates  plot  files  for  plot  types 
1  and  2.  PLTFI1.F0R  prompts  the  operator  for  the  scenarios  and  resources 
which  are  to  be  plotted,  using  the  five  character  scenario  code  previously 
defined.  PLTFI1.F0R  then  searches  the  archive  file  PSP1.ARC  for  the  resource 
and  scenarios  requested.  When  the  data  to  be  plotted  is  found,  it  is  written 
to  a  file,  whose  name  is  specified  by  the  operator,  in  a  plottable  ASCII 
format. 

Note:  When  generating  plot  type  4  the  first  scenario  entered  is  assumed  to 
be  the  reference  scenario. 
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1  !.  1  i  ;  0  ' ,  it  FORMATS  '  Hi'JUi  M A ,"}  Y  OEM AR I OS  TO  EE  p 

;  1 0  i  i  R  £  i .  0  S  5 . 6  0  ,  £  R  R  ~  4  0  ••  N  P  i._  0  T 

)  1  3  0  0  6  is  F 0 R  MAT  (  £  1  ) 

U  -V  0  0  DO  i  00  l~  l  .  i-JPLOT 


‘  'V;  H  0 

7  n 

TYPE  80  t 

;  6  •)  0 

p  n 

FORMAT ( "  TYPE  6  ■ 

■ ri  A P  h«.,  ;  p. 

7  ODE  F 

■  7  0  0 

P£hD( 3 , ? 0 v  ERR- 70 

.  CODES  C 

; 

1  3  n  o 

?o 

FORMAT v A* > 

-  ?  0  0 

1.  0  0 

CONTINUE 

2  ■  J  .•  Q 

.1  l  0 

7  y  p  i:”  P  0 

2 0  0 

i  ;•  n 

FORMAT  v  *  TYPE  THE 

:vnnE  C ;r 

THE  ruOT 

220  0 

R  £Ai>  (  5  9  0  £R  -■■■  i.  i  f;  ,• 

CUTF ::: 

230  0 

OPi:.N‘  UNI  F-3;  ACLE 

•=  EEOIN  ' 

FILE  ■■=  ■'  P 

2  >  0  0 

0  p  F.  N  (  U  N I T  a  4  A  C  C  E  3  3 

■'  8  E  0  0  U  T 

' , FILE : 0 

230  0 

WR  ITE  4 , 30  )  NAME 

2 0  0 

WR ITE< 4 .60 )  HP LOT 

27  1)  0 

IF  LOT- t 

28  0  0 

l  3  0 

RE ,40‘  3  90  .END  999° 

■  NCODE 

2'  0  0 

T  F  <  0  ■  0  ;>  E  (  £  P  L  0  7  )  ,  E  Q 

.  slC  DDE  > 

GOTO  140 

3  o  o  j 

COT a  130 

o  o 

•.  v  i; 

R E  Yi  0  ■’  3  i  5 0  4  MSG 

3  2  0  0 

i  3  0 

FORMA f ■ 20 A 4 ) 

33  0  0 

READ (3 ,  I  60  >  X  ,  NR 3 

4  o  o 

1  6  0 

F  OR  M A  I  (  1 3 > 4X  .  1 2  ) 

330  0 

READ < 3  ,  l 70 )  < NNAME 

(  I  >  >  I“l 

MRS) 

j  3  a  0  0  i  7  0  F  0  R  M  A  T  <  2  0  <  A  4 , 3  X  >  ) 

0  3700  INDEXED 

122 0  0  DO  180  J= 1 ; NR S 

1 3  9  0  0  I F  N  A  H  £  ,  E  Q  .  M  N  A  M  E  <  J )  )  I N  D  E  X  ~  J 

i)4i)l)0  [SO  CONI'  I  SUE 


4  l  0  0 

T  1“ 

(  INDE 

:< .  eq  . 

o  > 

GOTO  9990 

; ;■( 

;)  ;»•  ^ 

AD  •:  3 , 1. 

35)  N 

PER 

,  NKEPT 

■'*3  1  ('! 

•  33  FQi 

p  m  a  r (- 1 

6 . 3  X  , 

16) 

i .)  i) 

UR 

T.  TE  ■:  4  . 

1 35 ) NPER 

, NK EP T 

-TOO  0 

UR 

ITE  <  4  . 

150  ) 

MSG 

r  i’O  '  1  *  f 

1  ?  0  T 

*«  ,23 

READ- 3 

..  20  0  > 

S  c 

L  T  !..  ,  J  .;  .  i..  ~  1  ;  3  > 

AS";  (j  ■) 

1 no 

NT  I  HUE 

<1 1; 

20  0  EG 

F  HA  r  (  2 

0  v  2X  , 

E 1 1 

.  4  >  i 

3  '"i  0  0 

UR 

SITE  1  4  , 

21  0 ) 

(CL 

T  < INDEX . J >  ,  T - 1  ,23) 

3  0  i 

71  0  1-  i 

’  \  f  *. 

R  ivi  (A  T  <  2 

3  (  2  X  > 

HI  i 

,  4  )  ■ 

3  3  0  0 

....  /...  u  J 

REiYD  <  3 

-  ,  St 

,  20  0  > 

<  S 

I  MU  1.3  J  ,  L;  L-l  ,  Ni<  EP 

'  .  JOO 

22  0  CIO 

NT  T.  i'-'JE 

3  3  0  0 

UR 

ITE  1  1 

230  ) 

'■  S  i. 

i” l.J «...  (  INDEX  L  )  ,  L  ■■  1  .  N 

36  0  0 

2  30  FO 

rmiO  r  i ). 

o  o  <  ;?x 

;.  E  1 

1.  ,  4  )  ) 

370  0 

SO  IP 

LOT  I P 

...  0  T  +  1 

330  0 

R  E 

/i  I N  0  3 

3  ?  0  ) 

I  • ' 

■:  ipi.o 

r  .u: . 

i)  PL 

OT  cr.TO  1 20  A  Q 

a.Q 
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•  >  :•  -  , 


‘ii’";.  •  ij  /  ,  r;  SC  '  2 > 


G  •  :•  .  I  •  K . 


i 

FOPm, ,  r  •  V  .'  IE 

rC  SE  •' 

'•  0 

RvL  iO  s  ":  i;  ,  I;.:  p  ;  0  • 

FORMA  r  ■:  A 4  :• 

NAME 

-i  a 

TYPE  SO 

s  >) 

FORMAT.  '  HOW  ni.vi 

SEiiAR  IO 

P £AD  ‘  :i  .  6  ER  R  =4  0  • 

np  :..Q  r 

.  !  0 

6  0 

FORMAT. I i ■ 

'.  -V  :.i 

DO  I  TO  I  (  NRLCT 

i  -■)  u 

7  4 

TYPE  30 , I 

l  0  0 

SO 

FORMA  T.'  '  TYPE  6 

CHAR AC  1 

■  '  "0  !  1 

READ ( 5, 90, ERR -70 

>  CODE-: 

:  '?  i\l) 

9  0 

FORMAT  <  A.-> ) 

\  ?  )  ! ) 

S  0  o 

CONTINUE 

.  •'! 

!.  ,  0 

TYPE  120 

1  I 

■j  :j 

FORMAT-:'  TYPE  THE 

NAME  OF 

•  ,  ^  1  1  .) 

READ  ( 5  ,  ?(•  ,  ERR  =  HO  > 

OUTF  It... 

..  ^  J 

O  P  E  N  <  U  N I T  ™  3  ,  A  C  C  E  3  3 

-•  SEQIN 

: 

OP EN < UN I 7-4 , ACCESS 

-  '  SEQOU 

".  «  l  «.l 

WR  ITE  .'  4 , 30  ">  NAME 

:  i’.  ■) 

WRITE! 4 ,60 )  NPLOT 

•;  •  {) 

I PLOT  - 1 

'  >  0. 0 

1  3  0 

R  £  A  D  <  3 , 9  0  ,  E  N  D  *  9 9  9  9 

>  r>  CODE 

V>»J0 

IF  v CODE < I PLOT; , EQ 

.  NCODE ) 

5  ■)  <)  ) 

GOTO  >30 

>  1  :)  ) 

1  4  0 

READ ( 3 . 150)  MSG 

i  o  o 

l  •  4  0  0 

■ 'jar. 

i ;')  i)  Ij 

.)  ;,7!)o 
•L.1G0  0 

'  ?  a  n 

■  I  M) no 

•') !  i  n  o 
.  irr.no 

•  4.400 

o  jo  a 

i  ■  ’ o  o 

.  4.>n  n 
. -U>  i'i 
\  3  i'i  o 

i  [it) 

!  5  0  0  (1 
'■<" ;  ■’■) 
-  ">  z  • !  o 

' 1 ;n  i  i 
1 '  ’3  1 1 i  0 
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